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Abstract 

A model to study two-proton emission processes induced by electron scattering is developed. 
The process is induced by one-body electromagnetic operators acting together with short-range 
correlations, and by two-body A currents. The model includes all the diagrams containing a single 
correlation function. A test of the sensitivity of the model to the various theoretical inputs is done. 
An investigation of the relevance of the A currents is done by changing the final state angular 
momentum, excitation energy and momentum transfer. The sensitivity of the cross section to the 
details of the correlation function is studied by using realistic and schematic correlations. Results 
for 12 C, 16 and 40 Ca nuclei are presented. 



PACS number(s): 21.10.Ft, 21.60.-n 



1 INTRODUCTION 



This article belongs to a series dedicated to the study of the effects of the short range correlations 
(SRC) in processes induced by electromagnetic probes on atomic nuclei 

The SRC are produced by the strong repulsion of the bare nucleon-nucleon potential at internucleon 
distances smaller than ~ 0.5 fm. Every nuclear structure calculation that uses bare nucleon-nucleon 
potentials requires the presence of SRC, however, in medium- heavy nuclei, phenomena which can be 
unambiguously attributed to them have not yet been identified. 

A well known effect of the SRC is the depletion of the occupation probability of the quasi-hole 
states. Unfortunately, an analogous effect is also produced by the coupling of the single particle wave 
function with the low-lying phonons generated by collective nuclear vibrations. The inclusion of both 
effects seems to be necessary to account for the empirical occupation numbers ||. 

A more striking, and even better known, effect of the SRC is the various order of magnitude 
enhancement of the nucleon momentum distribution n(k) at high momentum values Q. Unfortunately, 
this quantity is not directly observable. It is therefore necessary to find measurable quantities related 
to it, hoping that they are sufficiently sensitive to n(fe) to allow for an unambiguous identification of 
the SRC effects. 

In these last years we have developed a model to describe the responses of finite nuclear systems 
within the framework of the correlated basis function theory. The model is inspired to the nuclear 
matter works of Refs. jj], ||]. With respect to these calculations, in our model the cluster expansion 
is truncated and only those terms containing a single correlation function are considered. The set of 
diagrams taken into account conserves the proper normalization of the many-body wave function. 

We tested the validity of this treatment by comparing the nuclear matter charge responses calcu- 
lated with our model with those obtained by the full expansion ||. The excellent agreement between 
these two results gave us confidence in extending the application of the model to other cases. Since 
our description of the nuclear excitations does not treat properly collective states, we have selected 
situations dominated by single particle dynamics. 

With our model we have studied the electromagnetic form factors of discrete states with large 
angular momentum [10]. These states are dominated by a single, or at most a few, particle-hole 
excitation. We found that the SRC produce small effects and are unable to account for the quenching 
required to reproduce the data. Electron scattering inclusive reactions in the quasi-elastic region have 
been investigated in Q. In this case the effect of the correlations is much smaller than that produced by 
the final state interaction, and therefore difficult to disentangle. We have also investigated phenomena 
where a single nucleon is emitted by the interaction with a virtual [^] or real || photon. The study of 
(e,e'p) reactions show a very small sensitivity to the SRC. The dominant effect beyond the independent 
particle model is produced by the final state interaction. Nucleon emission induced by real photons 
shows a larger sensitivity to SRC, especially for excitation energies around 200 MeV and large nucleon 
emission angles Q. Unfortunately, in this case, the meson exchange currents produce effects even 
larger than those induced by the SRC. 

The difficulty in identifying clear SRC effects in our previous studies was due to the presence 
of the large contribution of the uncorrelated one-body responses that dominates the cross sections. 
A possibility of eliminating the one-body responses is offered by two-nucleon emission experiments. 
Only two-body transition operators can induce the emission of two nucleons. These operators can 
effectively be constructed by one-body operators acting on a correlated nucleonic pair, or by meson 
exchange currents. If the two emitted nucleons are of the same type, two protons for example, the 
meson exchange currents contributions produced by the exchange of a single charged mesons do not 
contribute. These facts make two-proton emission experiments an ideal tool to investigate SRC. 

From the theoretical point of view, two-nucleon emission processes have been systematically studied 
in these years by the Pavia and Gent groups. The approach of the Pavia group [11, [l^] is based on the 
two-nucleon spectral function, and, from the theoretical point of view, it is a straightforward extension 
of the well tested approach used to describe single nucleon emission processes. Recently two-nucleon 
spectral functions calculated with microscopic theories [13|, 14] have been used. The Gent group [15j- 
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[p^tl , calculates transition amplitudes produced by effective two-body operators that are formed by 
connecting the traditional one-body electromagnetic operators to two-nucleon correlation function. 
The experimental situation is quite promising. After testing the feasibility of the experiment 



|Lq , angular correlations of the cross sections have been measured p0| , 21]. Other experiments of 
this kind have been planned. 

In this article we present the results obtained by applying our model to the study of two-proton 
emission induced by electron scattering. A comparison of our results with the available data is out 
the scope of the present article. Our aim is to investigate the sensitivity of the (e,e'2p) cross section 
to the details of the SRC. To this purpose we have studied the influence on the cross section of the 
various theoretical inputs required by the calculation. 

A summary of the basic formulae used to describe (e,e'2p) processes and electromagnetic currents 
is given in Section ||. In Section [3| we describe our nuclear model and in Section || we discuss the results 
we have obtained in 12 C, 16 and 40 Ca target nuclei. In the last section we draw our conclusions. 



2 THE CROSS SECTION 

A detailed derivation of the cross section for double coincidence experiments can be found in ]p2]] . 
In this section we briefly recall the expressions used in our calculations. We work in natural units 
(% = c = l,e 2 = 1/137.04) and employ the conventions of Bjorken and Drell |23|]. The initial and 
final electron four- vectors are respectively k = (e,k) and k! = (e',k'), and q = (u, q) = k — k' is the 
four-momentum transfer. The four-momenta of the emitted nucleons are indicated with p\ = (ei,pi) 
and p 2 = (e 2 ,P2)- 




Figure 1: Reference system used in our calculations. 

The reference system we adopt is shown in Fig. [l| The scattering plane is defined by the electron 
vectors k and k', 9 is the angle of the scattered electron and the quantization axis is taken along 
the direction of q. The directions of the emitted nucleons are determined by the angles ~f Pi and 4> Pi , 
i = 1,2. To simplify the formalism our calculations have been done considering two different angles 
ip Pi and 9 Pi . The relation to the angles described in the figure is 



and 

Opi -- 
9p2 '- 




b n if < <\> Pi < 7r 
b Pi - tt if 7r < <p Pi < 2ir 



Tpi if < 7 Pl < 7r 

if tt < 7 Pl < 2ir , 

if < 7 P2 < tt 
7p 2 if tt < 7p 2 < 2vr • 
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Note that the 7 P2 angle in Fig. |] has been defined positive in the half plane containing the scattered 
electron. The convention for 7 pi is opposite. 

We have derived the cross section expression by considering that the electron wave functions are 
plane wave solutions of the Dirac equation. We suppose that only one photon is exchanged between 
the electron and the nucleus and we neglect all the terms depending on the electron rest mass. With 
these approximations we obtain 

d 8 ^ K r / 

CMott /rec [VlWl + V t W t + V t lW t l + V U W U ) , (1) 



de'df^edeidf^df^ (2ir) 
where we have indicated with OMott the Mott cross section 



e 2 cos 



(0/2) \ 



2e sin 2 (9/2) J 



CMott : 

and with f Tec the recoil factor 

m lr»i I In I 

(3) 



r 1 =i+ m 

J rec x ' 



M A -2 



1 , l pi l a l q l a 
1 + i r COS P12 — i — r cos "p2 

|F»2 1 |P2| 



where #12 is the angle between the momenta of the two emitted nucleons, pi and P2, and m and Ma-2 
are the nucleon and rest nucleus masses, respectively. The factor K is 

if = m 2 | P i||p 2 |, (4) 

because we used non relativistic kinematics to describe the motion of the two emitted nucleons. The 
expression ([j]) has been obtained by integrating on the energy e 2 of one of the emitted particles and 
using the energy conservation. 

In the plane wave approximation the leptonic and hadronic vertexes can be factorized. In our 
expression (|l|) the terms v come from the leptonic vertex and depend only from kinematic variables 

vt = tan (6) 

VU = ( tan 2-#J ' (7) 

vtt = -4- (8) 
2q 2 

The information about the nuclear structure is included in the w factors. Because of the current 
conservation only three components of the four-vector current are independent. We choose the charge 
p(q) and the two transverse components in spherical coordinates 

J± = T^(J x ±iJ y ) ■ (9) 

In analogy to the (e,e'p) case ^ the w factors can be expressed as 

wt = <*i|pt(q)|^ / )($ / |p( q )|* i ) j (10) 
w t = (^|ji(q)|* / )(%|J_(q)|^) + (*i|4(q)|* / )(* / |J + (q)|* i }, (11) 
wu = 2ite((* i |pt(q)|* / )(* / |j_(q)|* i )_(* i |pt( q )|^ / )($ / |j + ( q )|^r i )) , (12) 

w tt = 2Re ((*<l4(q)|*/)(*/| J-(q)|*i)) , (13) 
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where we have indicated with \^fj) and |^/) the initial and final states of the full hadronic system. 
We do not consider the polarization of the emitted nucleons, therefore, in the previous equations, a 
sum on the third components of the spin of the emitted particles and of the angular momentum of 
the residual nucleus is understood. 

The one-body electromagnetic operators we have used are produced by the charge operator 



A 



1 + T 3 (fe) 



p(r) = £^^(r-r fc ), (14) 
k=i z 

and by the magnetization current operator 

^ M W - £ ^ (V 1 ^ + V x *(r - , 4 ) «(k) . (15) 

In the previous equations indicates the rest mass of fc-th nucleon, n P and the anomalous 
magnetic moment of the proton and the neutron, respectively, cr(k) the Pauli spin matrix of the fc-th 
nucleon and T^(k) = 1 (— 1) if the fc-th nucleon is a proton (neutron). In our calculations the nucleonic 
internal structure has been considered by folding the point-like responses with the electromagnetic 



nucleon form factors of Ref. [24|. To simplify the calculations we did not include the convection 



current because we know that its contribution is small in the quasi-elastic region 



7C 





7C 



A 



B 



Figure 2: Meson exchange currents diagrams considered in our calculations. The wiggly lines represent the 
virtual photon exchanged with the electron, the full thin lines the nucleons. Since we describe the emission of 
two protons, the exchanged pion is chargeless. 

Since we have considered only the emission of two like particles (two protons), the two-body 
currents induced by the exchange of charged mesons do not contribute. The main two-body current 
contribution arises from diagrams like those of Fig. Q where the exchanged it meson is chargeless. 
These diagrams correspond to the A-isobar currents for which we have used the following expression 



J (r) = 9m3 2^ r 3 {l)X7 x 5(r -r k ) 



k<l=l 



i2G A a(l) x V - 4G^V [or(fc) • VH(r - r/,ej)] 
+{k < — ► I), 



(16) 



where the factors f n NA, f-wNN and / 7 ata are the coupling constants related to the processes indicated 
by the lower labels, and m n is the pion mass. The Fourier transform of the above expression has been 
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multiplied by the A electromagnetic form factor that we have considered in its dipole form 



*A(|qM 



1 -2 



q- 



A 2 



(17) 



with A=855 MeV when u and q are expressed in MeV. The two coupling constants and are 
related to the two diagrams shown in Fig. ^ and are defined as 



G% = Ga ± Gb 



Following the indications of Refs. [25, 26], we have taken the coupling 

-l 

(J A = rriA 

where 



(18) 



(19) 



s = 2uj(m n + e h ) + (m n + e h ) 2 , (20) 
and e/j is energy of the hole single particle level. The expression of the imaginary part is taken from 



Ref. |7), 



Ta(V^) 



Unn (mA - m n ) (s - m 2 ) 3/2 



(21) 



3 ml y/s 
On the other hand, 

G B = (toa - mn)- 1 . (22) 

In the previous equations ttt-a = 1232 MeV and m n = 938.9 MeV are the A and nucleon masses, 
respectively. Finally, with the function T~C(r) we indicated the Fourier transform of the dynamical pion 
propagator 

d 3 fc F 7rN (A;,q)e ik - (r ~ r ' ) 



H(r - r;,e ; ) 



(23) 



(2tt) 3 k 2 + ml - e\ ' 

where is the pion-nucleon form factor and e\ is the energy of the exchanged pion obtained as 
the difference between the energies of the final and initial states of the Z-th. nucleon. To simplify our 
calculations we set the pion-nucleon form factor F 7T -^(k, e) = 1. We have verified that it is a very good 
assumption in the energy and momentum region we are interested to explore [EH]. 



3 THE NUCLEAR MODEL 

In the previous section we have presented the expressions of the cross sections and of the electromag- 
netic operators used in our calculations. Now we have to specify how we describe the initial and final 
nuclear states to be used to calculate the transition matrix elements in Eqs. (JlOj) - (|l3|) . To simplify the 
presentation we will be first concerned with the description of the final state which, asymptotically, 
has two nucleons in the continuum and two holes. In a following step we shall discuss how we describe 
the correlated many-body states. Our numerical calculations have been restricted to the case of the 
emission of two protons, but the formalism is general enough to deal with the emission of two generic 
nucleons. 

We describe the single particle wave function of a nucleon moving in the continuum with momentum 
p and with third components of spin and isospin s and t, respectively, as 

. oo Z -, 

Hp-*)XsXt = 7^ E E E {^) l RUpr)Y^{%)^-s\jm)Y^^ r ) Xt . (24) 
IPI l=o fi=-l j=l±l/2 

In the above equation we have done the multipole decomposition of the wave function and we have 
indicated with (lfj,^s\jm) the Clebsch-Gordan coefficients, with the spherical harmonics and with 

^(«) = E<*/Vb' m > y lA«( n )x- = l4i™>> ( 25 ) 
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the spin spherical harmonics [29|. The symbols \s an d Xt indicate, respectively, the spin and isospin 
parts of the wave function and R\ j{pr) its radial part, obtained by solving the Schrodinger equation 
with a spherical mean field potential. 

The many-body final state is characterized by the momenta of the emitted nucleons pi and P2, 
their spin and isospin components s±, S2, ii and t 2 , the total angular momentum of the residual nucleus 
J/i and its z-axis component M^, and the orbital and total angular momenta, I fa, jfa, I fa and jfa, of 
the single particle states where the nucleons are coming from 

1^(0,0) ) . (26) 

In the above equation a+, creates a particle with momentum |p| and spin third component s (the 
isospin indexes have been dropped to simplify the writing), and we have defined the annihilation 
operator as 

ai, J ,m = (-iy +m ai, j ,-m , (27) 

in order to deal with irreducible spherical tensors. The symbol (8> indicates the tensor product between 
two of such tensors 29]. 

The nuclear ground state is represented by |^(0, 0) ) defined as 

W ' 0)) -(*(0,0)|F+F|*(0,0))' (28) 

where F is the correlation function and |3>(0, 0)) is the Slater determinant describing a mean field 
ground state. We work with doubly magic nuclei, therefore the two zeros on the ground state symbol 
indicate the angular momentum value and its projection on the quantization axis. 

By considering the multipole expansion of the single particle continuum wave function we can 
write 

|*f(pi,Sl,tl,P2, s 2 ,t 2 ; Jh,M h ;lfa,jfa,l h2 ,jfa) ) 

-i£n£i E E E E 

'PlMpi 3v\ m V\ h 2 ^P2 JP2 m P2 



2 S l \jpi m Pl ) (lp2 \^P2 2 
(Jpi j /^pi ) Jp! > TTlpi j ^2 ' A%>2 ' JP2 ' m P2 > I hi 1 3 hi ! ^/l2 ' J^2 ! "^ft) ^h)) ■ (29) 



Yl pi m P1 (Opi) ^/ P2 m P2 (^2) (^PiMpi 9 s l|ipi m Pi) {IpzfJ'Pz r > s 2\jp2 m P2 



The many-body final state on the right hand side of the above equation is described as 

I VPl ) Mpi 5 Jpi ) m pi i ^P2 ' A*P2 > 3P2 J m P2 ) ^1 ' Jh\ ) ^2 ' ^2 ) M-h)) 

= E (-l)^ +mfc i + ^ +mk » Oh! -m^lAMfc) ^ r(30) 

The Slater determinant |$f) is built on the ground state Slater determinant |3>(0, 0)) by substituting 
the hole wave functions h\ and h 2 with the particle wave functions p\ and p 2 . In this specific case, 
the various single particle wave functions are not coupled to a total angular momentum, therefore the 
quantum numbers characterizing |3>f) are the orbital, I, and total, j, angular momenta of the single 
particle levels, their z-axis projections, [i and m, respectively, and the third components of the isospin. 

The evaluation of the w functions in Eqs. (|To|)-(13) requires the calculation of transition matrix 
elements of the type 

(¥ f |o„(q)|¥(o,o)) = (4 ^ 2 E E E E E ilpi+lp2 (-i) jh i +mh i +jh 2 +mh i 

' Pl ' ' P2 ' ipx/%1 3vi m Vl «p 2 Mp 2 ip2 m P2 m hl m h2 

{Ipi Mpi 2 S l Ijpi^pi ) (^P2^P2 2 S 2 b'p2 m 'P2 ) ~ ^2 ~~ m h 2 \JhM h ) 

y* P1 m pl (Shi) YC P2mp2 (n p2 ) e[0,(q) : Pi, P2 ,h 1 ,h 2 ] , (31) 
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where 



0„(q) 



<Tre 



-iqr 



O v (r), 



indicates a one-body electromagnetic operator and we have defined 

(cD f |F+O^(q)F|ci>(0,0)) 



£[C*r,(q) :px,P2,hi,h 2 ] 



($ f |F+F|$ f )2($(0,0)|F+F|$(0,0)) 



(<f> f \F+F\<S> { ) 



($f\F+F\<f>f) 



($(0,0)|F+F|$(0,0)> 



In our calculations we have used purely scalar correlations defined as: 



(32) 



(33) 



F(l,2,...A)=Uf(, 



i<j 



where r« = |rj — r ^ | is the distance between the positions of the particles i and j. 



(34) 



The two factors in Eq. ( p3[ ) are separately evaluated by expanding numerator and denominator in 
powers of the two-body short-range correlation function /. The presence of the denominator and the 
energy conservation eliminate the unlinked diagrams. 

Since in our calculations the correlation functions are purely scalar they commute with the operator 
Ojj(q) therefore: 



£[Or/(q) ■Pl,P2,hl,h 2 ] 



($;p 1 ,p 2 ,h 1 ,h 2 \O v (c L )F 2 \<I>;00,+l) L 

A 

{$;pi, P2 , h lt h 2 \O v (q) n(l + M|$(0,0)) L , 

i<3 



(35) 



where we have used the function h+j = f 2 (rij) — 1 and the subindex L indicates that only the linked 
diagrams are considered. 

The approximation of our model consists in retaining only those terms where the function 
appears only once 



£[O v {q) : pi,p 2 ,h 1 ,h 2 ] 



^[O v (q) : pi,p 2 ,hi,h 2 ] 

(*;pi,pz, hi,h 2 \ 0,(q) + h v) °)>£ 



i<j 



{&,Pl,P2, h 2 \ 0„(q) h (ri,j) MO, 0)) L 

+ ($;p 1 ,p 2 ,h 1 ,h 2 \O r ,(q) J2 h Kj)\$(0,0))L- 

Ki<j 



(36) 



(37) 



This result has been obtained using a procedure analogous to that adopted in Ref. |30| for the 
evaluation of the density distribution, therefore the truncation of the expansion is done only after the 
elimination of the unlinked diagrams. The Meyer- like diagrams considered in our work are presented 
in Fig. ||[ It is evident that, in this case, the uncorrelated term, produced by the 1 in Eq. (pq) , does 
not contribute. All the terms necessary to calculate the cross section have been specified. The detailed 
calculation of the various matrix elements for the one-body current is presented in the Appendix A. 

The A currents allow for the two-nucleon emission already at the independent particle model 
level. We consider only this contribution, therefore SRC and A currents are not directly coupled in 
our model. They influence each other only via the interference terms. 



4 SPECIFIC APPLICATIONS 

The interest in studying two-nucleon emission cross sections, is related to the possibility of obtaining 
information on the SRC. In order to identify correlation effects it is necessary to keep under control 
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(2.1) 



(3.1) 



(3.5) 



(2.2) 



(2.3) 



(2.4) 




1 P2 



+ 





2 Pl\ A 2 Si P 



+ 




(3.2) 



(3.3) 



(3.4) 



j 2 P ,\^2 j>2 P2 \j fl j> p 2 y 2 



+ 



v' + V " *v 



(3.6) 



(3.7) 



(3.8) 



Figure 3: Meyer-like diagrams representing the one-body currents matrix elements included in our calculations. 
The dashed lines represent the function h(r) = f 2 (r) — 1, where f(r) is the correlation function. The black 
squares are the points where the electromagnetic one-body operator O^(r) is acting. The oriented full lines 
represent single particle wave functions and one has to understand a sum on all the hole states labeled with i. 



all the other variables which can affect the cross section. For this reason, before studying the influence 
of the SRC we shall discuss the sensitivity of the cross section to the various theoretical inputs of our 
model. The results of this investigation will be presented only for the reaction 16 0(e,e'2p) 14 C, but 
we have done calculations also for 12 C and 40 Ca nuclei. The two-hole composition of the various final 
states we have considered is given in Table [l]. 







16 


4U Ca 


of 


(lp3/2)- a 


(lpl/2)- 2 


(ld3/2)- 2 


0+ 




(lp3/2)- 2 


(2sl/2)- 2 


1+ 




(lpl/2)- 1 (lpS/2)- 1 


(ldS/2)- 1 (2S1/2)- 1 




(lp3/2)- 2 


(lpl/2)- 1 (lpS/2)- 1 


(ld3/2)- 2 


2+ 




(lp3/2)- 2 


(ldS/2)- 1 (2S1/2)" 1 



Table 1: Two-hole compositions of the final states considered in our study. 
4.1 Kinematics and No-Recoil Approximation 

The definition of the kinematics of the reaction under investigation is far from being trivial. Energy 
and momentum conservation imply 



u = ei + e 2 - e hl - e h2 + 

q = Pl + P2 + Pr , 



Pr 



2M A - 



(38) 
(39) 
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where we have considered non relativistic kinematics and we have indicated with p r the recoil mo- 
mentum of the A — 2 residual nucleus. The other symbols have been already defined in the previous 
sections. 

In our calculations we fixed oj, q, and P2. By selecting the A — 2 nucleus final state also 
and eh 2 are fixed. The cross sections are presented as a function of 7 Pl . By solving the above 
system of equations we obtain the values of |p r | and |pi| for each r y pi . This means that in the 
calculation of the nuclear transition matrix elements, the radial integrals shown in Appendix A, 
should be calculated for each value of j pi . From the computational point of view, the situation 
improves enormously by assuming p r = 0, since in this case |pi| becomes independent from j pi . We 
call No-Recoil Approximation (NRA) this approximation, consisting in setting to zero the kinetic 
energy and the momentum of the residual A — 2 nucleus. 



16 O (e,e'2p) 14 C 




-180 -120 -60 60 120 180 
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Figure 4: Cross sections for the 16 0(e,e'2p) 14 C reaction in the standard kinematics (see Sect. h2). The 
two-hole composition of the final states is given in Table [|. The dashed lines have been obtained by properly 
considering the recoil energy of the 14 C nucleus. The full lines have been obtained within the NRA. The left 
and right panels show the same results in linear and logarithmic scale. 

In Fig. | we compare the results of the NRA (full lines) with those obtained by correctly treating 
the A — 2 nucleus recoil (dashed lines). The reaction is 16 0(e,e'2p) 14 C and the calculations have been 
done in what we shall define in Sect, [h^ as standard kinematics. The left and right panels show the 
same results in linear and logarithmic scale. This presentation is done to show that the two results 
are very similar in the peak region. The differences are relevant only off the maximum region, but 
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they can be appreciated only in logarithmic scale since in that region the cross section values fall by 
several order of magnitude with respect to the peak values. 

We have verified that, in the kinematics used in the calculations, the recoil factor in Eq. ([|) differs 
from 1 at most by a 4%. Therefore the source of the difference is related to the change of the single 
particle wave functions because of the change of ei = pf /2m. In the kinematics under investigation 
this energy can vary from 34.0 down to 22.6 MeV. It is interesting to notice that, in the region where 
the cross section has its maximum, or maxima, the results of the two calculations are rather similar. 
This feature can be more easily understood by considering the emitted nucleon wave functions as 



plane waves. In this case by using Eq. (|32D one can observe that the calculation of the matrix element 



(35) is related to the Fourier transform of the two-hole relative wave functions with respect to the 



variable p r = q — (pi + P2). As it has been pointed out in Refs. [12, 16, |17]], the maximum value of 
this transform occurs when the argument has its minimum value, possibly zero. For this reason the 
results of the two calculations shown in Fig. || are rather similar in the region of the maxima. 

The good performances of the NRA, at least in the region where the cross section shows its maxima, 
gave us confidence about its use in our calculations. All the results presented in the following sections 
have been obtained in NRA. 

4.2 The mean field 

In this section we present the results of our study on the influence of the mean field used to generate 
the single particle wave functions on the (e,e'2p) cross sections. The investigation has been conducted 
on 16 O by fixing the correlation function and changing the mean field. 

As in the case of our previous works on (e,e'p) and (7,p) reactions || ||, we used a real potential 
to generate the hole wave functions, and an optical potential to generate the particle wave functions. 
With this choice the single particle wave functions are not any more orthonormal. In our previous 
investigations [§, [| we have studied these non-orthogonality effects, and we reached the same conclu- 



sions obtained in Ref. [31] where the problem has been first studied: for the case under investigation 
and the kinematics we intend to explore, these effects are negligible. 

The parameters of the Woods-Saxon well used to generate the hole wave functions are those given 



in Ref. [32]. Associated to these single particle wave functions there are SRC obtained by minimizing 



the nuclear hamiltonian with the semi-realistic S3 interaction of Afnan and Tang 33]. For our test 



cases we used the gaussian form of the correlation function taken from 32]: 

f(r) = l-ae~ br \ (40) 

with a=0.7 and 6=2.2 fm -2 . The particle wave functions have been obtained by using the optical 
potential of Schwandt et al. p4j]. If not explicitly stated otherwise, these are the parameters used in 
our calculations. 

Since the number of variables is quite large we have chosen specific kinematics to perform our 
test calculations. These kinematics, we call them standard kinematics, consist in fixing the electron 
incoming energy e = 800 MeV and the excitation energy u = 100 MeV. Then we fixed £2 = 40 MeV 
and the emission angle 7 P2 = 60°. The energy of the other proton is obtained from Eq. (|38| ) by setting 
p r = 0, that is in the NRA, and the cross section are presented as a function of j pi , the emission 
angle. We always work in coplanar kinematics, i.e. 4> Pl and (p P2 are zero. We recall that we have taken 
the angle j pi to be positive on the opposite side of the scattered electron with respect to q, while 7 P2 
is taken positive on the same side of the scattered electron. We used this standard kinematics in 16 
as a reference calculation. 

A first test has been done on the sensitivity of the results to the changes of the hole wave functions. 
An example of the results we have obtained is shown in Fig. ||. The standard results are shown by the 
full lines. The dashed lines have been obtained with the hole wave functions generated by the Woods- 



Saxon potential whose parameters are given in Table V of Ref. [35|. These parameters have been used 
in a Fermi hypernetted chain calculation done with the Argonne V8' nucleon-nucleon interaction. In 
this case the spin orbit term of the mean-field potential has been set to zero. The dotted lines have 
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Figure 5: Cross sections calculated in standard kinematics with different mean field potentials generating the 
hole wave functions. The full lines show the standard results obtained by using the mean-field potential of Ref. 
|32). The dashed lines have been obtained with the mean field potential of Ref. J|^]. The dotted lines are 
results obtained with harmonic oscillator wave functions. 



been obtained by using harmonic oscillator single particle wave functions. The sensitivity of various 
results to these changes is limited: we found a maximum variation of 4%, except for the case of the 
1 + state. This state is, however, out of the systematics, since it is dominated by the A currents, as 
we shall discuss in more detail in Sect. 4.2. 

We should remark that all the mean-field potentials used in these test calculations reproduce the 
root mean square charge radii. Furthermore, the single particle energies have been kept constant in 
all the calculations. In our calculations, the hole single particle energies are used to determine the 
particle energies, see Eq. (|38[). The particle wave functions are calculated by solving the single particle 
Schrodinger equation with an energy dependent optical potential for a continuum wave. Changes in 
the energy of the particle could strongly modify the single particle hamiltonian, then the particle wave 
functions. 

In our standard kinematics calculations we used the hole energies obtained by solving the single 
particle Schrodinger equation with the ground state Woods-Saxon potential. In the cases of our interest 
these single particle energies are given in Table |2[ In this approach we do not consider the fact that a 
residual nucleon-nucleon interaction rearranges the mean field in the A — 2 nucleon system such that 
its energy is not any more the energy of the A nucleon system minus the two single particle energies. 
A way to take into account this fact is to use the experimental binding energies of the A — 2 nucleus, 
the 14 C in the specific case, and the energies of the first excited states which should be interpreted in 
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ld5/2 


2sl/2 


lpl/2 


lp3/2 










-18.0 


16Q 






-12.7 


-16.9 


40 Ca 


-8.7 


-19.3 







Table 2: Single particle energies, in MeV, for the states of interest in our calculations. These energies have 
been obtained using the Woods-Saxon parameterization given in |32j . 



terms of single particle excitations. This is the approach used by Giusti and Pacati [12] 
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Figure 6: Comparison between cross sections calculated in standard kinematics (full lines) and those obtained 
by increasing (dashed lines) and reducing (dotted lines) t\ by 2.5 MeV. 

From the pragmatic point of view the two methods of determining the hole single particle energies 
introduce an uncertainty of about 2.5 MeV. This uncertainty is reflected on the energies of the emitted 
particle and on their wave functions. The final uncertainties are shown in Fig. |6| where the results 
of the standard calculations are compared with those obtained by increasing and lowering ei by 2.5 
MeV. 

The larger variation of the cross section obtained in Fig. ^ with respect to those of Fig. |5| indicates 
the great sensitivity of the results to the particle wave functions. In Fig. ^ we compare the results 
of the standard calculation with those obtained by using plane waves for the particle wave functions 
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(dotted lines) and with those obtained by using the same real Woods-Saxon potential utilized to 
describe the hole states (dashed lines). It is evident the large difference between the various results 
and it appears clear the effect of the complex optical potential whose imaginary part quenches the 
cross section. 




l Pl [deg] 



Figure 7: Cross sections calculated in standard kinematics. The standard results (full lines) are compared 
with those obtained by using plane waves (dotted lines) and real Woods-Saxon single particle wave functions 
(dashed lines). 



The sensitivity of our calculations to the choice of the optical potential is shown in Fig. |8| where 
the cross sections have been calculated by using the optical potentials of Schwandt et al. |p3 (the 



standard calculations), Nadasen et al. [36| and Comfort and Karp [37|. Though the modifications of 
the cross section are less pronounced than in Fig. [5], they can become as large as 18%. 

We can conclude this section by saying that the various plausible choices of the mean-field param- 
eters can produce uncertainties of about 20% - 30% on the maxima of the cross sections. 



4.3 The two-body current 

We have already pointed out that the presence of SRC is not the only mechanism inducing the 
emission of two protons. Meson exchange currents of the type shown in Fig. [2] can also contribute to 
this process. For the purpose of studying SRC the presence of these two-body currents is disturbing, 
therefore, it is important to find kinematics where the A currents contributions are negligible with 
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Figure 8: Cross sections calculated in standard kinematics. The standard results (full lines) are compared with 
those obtained with the optical potential of Comfort and Karp |37J (dotted lines) and that of Nadasen et al. 
|36| (dashed lines). 



respect to the SRC effects. 

In this section we discuss the sensitivity of our results to the presence of the A currents. There 
are various aspects of the description of electromagnetically induced A excitation that are still largely 
debated [38, 38]. It is not our aim to enter in this discussion. We have considered a A current 
model commonly used in the literature, we have applied it to two-nucleon emission processes, and we 
have analyzed the sensitivity of our results to the changes of the parameters of the model and of the 
kinematics. 

We present in Fig. ^ some specific result obtained in standard kinematics for uj = 100 MeV and 
u = 200 MeV with the aim of discussing two different features regarding the A currents. The first 
one is related to the sensitivity of the results to the various coupling constants defined in Eq. (1C). 
The second one concerns the relative importance of SRC and A currents with respect to the change 
of the excitation energy. 

Out of the three constants f n NN, fwNA and / 7 jva> only the first one is rather well defined by the 
experimental data. Its commonly accepted value is f% NN /4tt=0.079. The values of the other two 
constants are not so precisely known. In our standard calculations we used the values taken from [flc([. 
In Figs. H we have also used the values adopted in Refs. 12] and [17] (see Table ||). 

We have also tested the validity of the static approximation in the description of the A propagator. 
This approximation consists essentially in setting Ga = Gb = l/(?nA — m n) in Eq. (|^). This is 
the approximation we have used in our previous works []l], |3|, [|] . The dotted lines of Fig. || show the 
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Figure 9: Cross sections for the 16 0(e,e'2p) 14 C process calculated using different values of the A current 
coupling constants (see Table ||). The full lines have been obtained with the parameters of Ref. and these 
are the values used in our standard calculations. The dashed and dashed dotted lines have been obtained with 
the parameters of Refs. |l2) and |l7[] respectively. Dotted curves correspond to the standard calculations but 
using the static propagator for the A currents. 





AMA GIU RYC 


/7JVA 
fnNA 


0.299 0.373 0.120 
1.69 2.15 2.15 



Table 3: Values of the parameters used in Eq. (|l6|). The AMA, GIU and RYC values are from Refs. g(J, 
and respectively. 



15 



results obtained with the static approximation. 

We observe that at u = 100 MeV the Of state shows scarce sensitivity to the changes of the values 
of the A currents coupling constants. More sensitive is the 2f state, but certainly the most sensitive 
one is the 1 + . These different degrees of sensitivity are connected to the relative importance of the 
A current contribution with respect to the SRC in the various excited states. As already pointed out 
in Refs. |l^] and [17], the 1 + state is strongly dominated by the A currents, therefore any change in 
the coupling constants produces large modifications of the cross section. The small sensitivity of the 
Of state is due to the fact that its excitation is mainly due to the longitudinal response (see Fig. 11), 
while the A acts on the transverse response only. The A contribution is not negligible for the 2f state 
and it becomes dominant for the 1 + . 

The results of the static approximation have been obtained by using the constants of Ref. [flC|]. 
This approximation is valid when its results are closed to those shown by the full lines. At lo=100 
MeV this is the case for all the states we have investigated but the 1 + . We do not show the results 
for the Of and 2f states that present analogous features of 0^ and 2f results. 

The situation changes for the excitation energy lo=200 MeV. Now we are in the resonant region 
of the A and here all the states are strongly sensitive to the different values of the coupling constants. 
The transverse response dominates also in the 0^ state. The results of the static approximations are 
quite far from the full lines, clearly showing the inadequacy of this approximation. It is interesting 
to notice that the static approximation produces results rather similar to those obtained with the 
coupling constants of Ref. [17]. 

We have investigated the energy and the momentum dependence induced by the SRC and by the A 
currents. The results of our study about the energy dependence are summarized in Fig. 1C, where in 
the left panels we show the cross sections calculated using only the one-body currents plus SRC, while 
in the right panels we show the results obtained with the A currents alone. The comparison between 
the left and right panels for the 1 + state, clearly show that the calculations with the A currents 
produce cross sections much larger than those obtained with one-body currents only. This confirms 
that the excitation of this state is dominated by the A currents dynamics. The same analysis requires 
a more careful discussion in the case of the 0^ state. At u = 100 and 150 MeV, the cross sections 
obtained with the one-body currents are larger than those calculated with the A currents alone by a 
factor 20-30. At oj = 200 MeV the difference between the two results is about five times smaller. This 
indicates that, even for the Of state, the contribution of the A currents becomes important when the 
excitation energy is in the A resonance region. For the 2f state the contributions of the one-body 
and of the A currents have the same order of magnitude. 

The SRC do not show specific energy dependence. We have verified that by calculating all the 
16 final states of Table |l[ In the example shown in Fig. |H] one can observe that in the 0^ state the 
cross sections at 200 MeV is smaller than at 150 MeV, while the size of the cross sections for these 
excitation energies is almost the same for the 2f state. The situation is different for the case of the 
A currents whose contribution grows with the energy since it enters in its resonance region. 

Some results of our investigation of the momentum transfer dependence are shown in Fig. 11. In 
this figure, instead of the cross sections, we present the longitudinal and transverse responses, Eqs. 
(10) and (11), to get rid of the trivial q dependent terms in the kinematics factors. We show the 
responses for the 0^ and 1 + final states. 

The Of state is dominated by the longitudinal response. The shape of the angular distribution of 
wi is not modified by the change of |q|, however its peak value is reduced with increasing |q|. The 
values of wt without A currents are so small to be negligible. The inclusion of the two-body currents 
enhances the transverse response. It has been shown in Ref. []14|] that, in this state, the emission of the 
two nucleons is dominated by a relative s-wave. The |q| dependences of longitudinal and transverse 
responses are quite different. At |q| = 300 and 400 MeV/c, the longitudinal responses are much larger 
than the transverse ones. The peak values of the responses become smaller with increasing |q|, but 
the decrease of u>i is steeper than that of wt- At |q|= 600 MeV/c we are in a situation where the size 
of the two responses is about the same. In this kinematics the presence of the A currents noticeably 
affect the full cross section. 
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Figure 10: Cross sections for the 16 0(e,e'2p) 14 C process calculated for three values of the excitation energy. 
Full lines w=100 MeV, dotted lines w=150 MeV, dashed lined w=200 MeV. The left panels show the results 
obtained by considering only the one-body currents plus the SRC. The right panels show the cross sections 
obtained with A currents only. 



The situation is different for the 1 + state. This state is dominated by the A currents, therefore 
the transverse response is the most important one. The emission of the two protons is ruled by a 



relative p-wave [14], which produces an angular distribution with two peaks. These two peaks are 



clearly visible when the kinematics verifies q ~ Pi + P2- If q S> Pi + P2 for any value of j pi , the two 
peaks merge as it is shown by the dashed and dashed dotted lines. The change of |q| does not modify 
sensitively the size of the peaks. 

To summarize the results of this section, we may say that, in order to minimize the effects of the A 
currents, one has to work at excitation energies far from the resonance peak, at relatively low values 
of the momentum transfer (|q| < 500 MeV/c), and furthermore one should not consider the 1 + state. 



4.4 The Short Range Correlation 

After having studied the sensitivity of our results to the theoretical inputs and to the presence of the 
A currents, we arrive now to the main point of our investigation: the sensitivity of the cross section 
to the details of the correlation function. It is obvious that the use of very different correlations 
would produce rather different results, since f(r) is the key ingredient of the calculation. The point 
is to investigate the sensitivity of the cross section to relatively small changes of realistic correlation 
functions. 

For our investigation we have considered the purely scalar correlation functions f(r) shown in Fig. 
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Figure 11: Longitudinal (wi) and transverse (wt) responses, Eqs. (Q and (|ll|), calculated in standard 
kinematics, for various values of the momentum transfer: dotted lines 300 MeV/c, full lines 400 MeV/c, dashed 
lines 500 MeV/c and dashed-dotted lines 600 MeV/c. The two central panels have been calculated without A 
currents contribution, while the curves of the lower panels include also this contribution. 



|l~2^ . The correlation function G, we have used up to now, and the S3 correlation function, are taken 
from Ref. J3^ ]. They have been fixed by minimizing the energy functional calculated with a nuclear 
hamiltonian containing the Afnan and Tang nucleon-nucleon interaction ]33[| . The two minimizations 
have been done by considering respectively a gaussian type correlation (G) and the Euler procedure 
described in that reference (S3). The correlation labeled as V8 in the figure, is the scalar part of the 
state dependent correlation used in Ref. 1 35 ] , where the hamiltonian was built by using the nucleon- 
nucleon Argonne V8' interaction plus the Urbana IX three-body force. The three correlation functions 
differ for few details. The 53 and V8 overshoot the asymptotic value of 1 in the region between 1 and 
2 fm. The V8 correlation function has a lower minimum than the other two. 



In the Figs. [L3, 14 and [15] we show the cross sections calculated in standard kinematics with the 



three correlation functions shown in Fig. 12 for three target nuclei: 12 C, 16 and 40 Ca. The list of 



the two-hole wave functions describing the final states of the A — 2 rest nuclei is given in Table |. 

Apart from details related to the specific nucleus and the final state considered, the general trend 
shown by these results is that the shapes of the angular distributions are not modified by the correla- 
tions. It is rather the size of the cross sections in their maxima which has changed. This indicates that 
the shapes of the angular distributions are ruled by the kinematics and by the angular momentum 
coupling between the hole states and the allowed partial waves where the two particles can be emitted. 
This has been widely discussed in [[lj] where a decomposition of the partial wave in terms of center 
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Figure 12: Correlation functions used in our calculations. The line labeled G indicates the gaussian correlation 
function of Ref. (3^] and that labeled 5*3 the Euler correlation function of the same reference. With V8 we 
show the scalar term of the state dependent correlation function of Ref. |}5| . 

of mass and relative motion of the emitted nucleon pair has been done. 

The cross sections obtained by using the V8 and S3 correlations are rather similar and they are 
roughly a factor two smaller than those obtained with the gaussian correlation, which we took as 
standard results. This effect is certainly larger than the theoretical uncertainties estimated in section 
fl.2j . The 1 + states are out of systematics because they are dominated by the A currents. 
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Figure 13: 12 C(e,e'2p) 10 B cross section in the standard kinematics for the emission of the two proton from the 
lp3/2 level of the target nucleus. The full, dotted and dashed lines, have been obtained by using respectively 
the G, S3 and V8 correlation functions of Fig. [12. 



To investigate the sources of these differences we have done calculations with rather schematic 
correlations which allowed us to switch on and off various characteristics of the correlation functions. 
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Figure 14: 16 0(e,e'2p) 14 C cross sections calculated in standard kinematics with the three correlation functions 
shown in Fig. 



To this purpose we used step function correlations of the type 

( A, < r < di 
f(r) = \ B, d 1 <r<d 2 (41) 
[ 1.0, d 2 < r 

We have changed the values of the A,B,d\ and d 2 parameters to produce four different correlations. 
The values of the parameters used in our calculations are given in Table |I[ 





A 


B 


di [fm] 


d 2 [fm] 


CSC1 


0.2 


1.0 


1.0 


1.0 


CSC2 


0.1 


1.0 


1.0 


1.0 


CSC3 


0.2 


1.1 


1.0 


1.3 


CSC4 


0.2 


1.0 


0.9 


0.9 



Table 4: Parameters of the schematic correlation functions as given into Eq. |4l]. 

The correlation function CSC1 is a square well of radius 1.0 fm and depth 0.2. The 16 0(e,e'2p) 14 C 
cross sections calculated in standard kinematics with this correlation function are shown in Fig. |l6| by 
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Figure 15: 40 Ca(e,e'2p) 38 Ar cross section calculated in the standard kinematics. The full, dotted and dashed 
lines, have been obtained by using respectively the G, 5*3 and V8 correlation functions of Fig. [13. 



the full lines. We shall discuss the differences produced by changing the correlations with respect to 
these results. The correlation function CSC2 has the same shape of CSC1 but a lower minimum. The 
effect is a small increase of the cross sections (see dashed curves). The CSC3 correlation function has 
a square shape of the type of that of CSC1 but, for 1.0 fm< r < 1.3 fm, it overshoots the asymptotic 
value of 1.0. This modification in the correlation function is enough to lower to half, roughly, the cross 
section maxima (dotted curves). An analogous effect is obtained by the CSC '4 which has the same 
shape as CSC1 but a smaller healing radius. The results obtained with this last correlation function 
are shown by the dashed-dotted curves. 

The behaviors just described can be understood by remembering that the quantity entering in 
the cross section calculation is not f(r) but rather h(r) = f 2 (r) — 1. The largest is the contribution 
of h(r) in Eq. fl37|), the largest are the cross sections. The CSC2 correlation has a larger overlap 
with the other integrated functions than CSC1, since its depth is deeper. The overshooting in CSC3 
generates a term in h(r) of opposite sign with respect to the rest of the function, therefore the total 
contribution to the integral becomes smaller. The same effect can be obtained by reducing the range 
of the correlation, as it is done in CSC '4. 



It is now easier to understand the results of Figs. 13J, |14| and [15. The three correlations used in 
the calculations reach their asymptotic values at about the same internucleon distance: 2 fm. The S3 
and V8 correlation functions overshoot this value in the region between 1 fm and 2 fm. This part has 
opposite effects as the part below the value of 1, as it has been showing the case of CSC3. Since these 
functions should be multiplied by a r 2 in the integral, the effect of the overshooting is larger than the 
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Figure 16: Cross section in standard kinematics calculated with the schematic correlations of Eq. (|4l|) using 
the parameters given Table f|. 



small differences in the minima at r = fm. 



It is remarkable that very schematic correlation functions as those used in Fig. 16 produce angular 



distributions with the same form as those of Fig. 14, obtained with realistic correlation functions. This 
confirms again that the shapes of the angular distributions are determined by the quantum numbers 
of the final states. 



We have tested the role of various diagrams considered in our model. In Fig. 17 the cross sections 
calculated in the standard kinematics and with the gaussian correlation function are presented. The 
dotted lines have been obtained by using the 2-point diagrams of Fig. ||. The dashed lines show the 
results we got when 2 and 3-point diagrams of Fig. |3| are used. The total results when also the A 
current contributions are added are presented by the full lines. 

The inclusion of the 3-point diagrams lowers the cross sections calculated with the 2-point diagrams 
only in all the cases we have considered. This reduction is not constant but it depends from j Pl , as it 
is clearly shown by the two 2 + states where the peaks around J P1 =0 are more quenched than those 
around 7 Pl =100°. The effect of the A currents does not have a regular behavior. The global quenching 
of the cross section shown in 0^ is followed by and enhancement in the 0^ state not shown in the 
figure. Also the two 2 + states show different behavior of the A currents. As already mentioned the 
1 + is dominated by the A currents even at 100 MeV. 

We have already mentioned that p r = q — pi — P2 is a crucial quantity in the two nucleon emission 
process. We have changed this quantity in our standard kinematics by modifying the j P2 angle. In 
Fig. 18 we present some of the results we have obtained using 7 P2 =30° (dotted lines), 7 P2 =60° (full 
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Figure 17: Two nucleon emission cross sections calculated in standard kinematics. The dotted lines have been 
obtained considering only the 2-point diagrams of Fig. ||. The dashed lines have been obtained by adding 
the 3-point diagrams and the full lines show the results when also the A currents are included. The gaussian 
correlation of Fig. [12 has been used. 



lines) and 7 P2 =90° (dashed lines). 

It is interesting to notice that the angular distributions of the + states do not change their form 
by changing ~f P2 , only the size of the peak is modified. The situation is rather different for the 2 + 
states. By changing 7 P2 there is a merging of the two peaks which are well separated at 7 P2 =60°. This 
situation is analogous to that discussed in the description of Fig. 11 



5 SUMMARY AND CONCLUSIONS 

We have studied electron induced two-nucleon emission processes by applying a model which has been 
already used to describe inclusive |l], [l^] and one nucleon emission processes 0, M. In this model, 
inspired to the correlated basis function theory, the full cluster expansion is truncated to consider all, 
and only, those diagrams containing a single correlation function h(r). 



In addition to the traditional two-point diagrams, evaluated also in other approaches [17|], we 
consider also three-point diagrams. These diagrams are necessary to conserve the proper normalization 
of the nuclear final state. The presence of these tree-point diagrams always reduces the two-point cross 
sections. From the quantitative point of view, in the kinematics we have investigated, we found relevant 
effects only for the 2 + final states, where we could detect differences of about 40%. In all the other 
cases the role of the three-point diagrams is much smaller. 
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Figure 18: Cross sections calculated in the standard kinematics but changing the value of the 7 P2 angle. Dotted 
lines 7 P2 =30°, full lines 7 P2 =60° and dashed lines 7 P2 =90°. 



We have shown two-proton emission cross sections as a function of the angle of one of the emitted 



protons. Because of the energy and momentum conservation, Eqs. (|38| ) and (39), this kinematic 
set up implies the change of the energy of the proton for every emission angle. We have used an 
approximation consisting in neglecting the recoil of the A — 2 residual nucleus, therefore the energy of 
the emitted nucleons is the same for every emission angle. This approximation works well in the region 
where the cross sections have their maxima. This indicates that these maxima appears when the recoil 
momentum of the residual nucleus has its minimum value, as was already discussed in [12, |l^, [17]]. 
Since the cross sections values drop by several order of magnitude in the regions outside the maximum, 
we used this no-recoil approximation in our study. 

We have investigated the sensitivity of our results to the possible changes of the mean-field input. 
We have modified the mean field parameters such as the description of the nuclear data different 
from (e,e'2p) reaction would not be compromised. For example, the hole wave functions have been 
modified but we always cared to reproduce the charge root mean square radii. We used different 
optical potentials taken from the literature but all of them reproduce elastic scattering proton-nucleus 
data in different energy regimes. We have estimated that plausible changes of the mean-field bases 
would imply a theoretical uncertainty of about 20% - 30% on the maxima of the (e,e'2p) cross section. 

The cross sections are most sensitive to the changes in the particle wave functions. Calculations 
done considering these wave functions to be plane waves, or by using purely real Woods-Saxon waves, 
show cross sections a factor two larger than those obtained with a complex optical potential. A 
consistent description of the one nucleon emission processes H |3| requires the use of complex optical 
potentials to take into account the excitation of emission channels different from the one considered. 
The effects of changing optical potential are certainly smaller than those produced by the use of real 
mean-field, but they are not irrelevant. In the kinematics we have studied we found changes as large 
as 18%. 

The presence of the two-body A currents interferes with the SRC in the two-nucleon emission 
mechanism. To minimize the effects of the A currents one has to find proper kinematics conditions. 
We found that the role of these two-body currents at excitation energy w=200 MeV is always relevant 
for every final state. This is reasonable since 200 MeV is well inside the A resonance region. 
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The situation changes for u>=100 MeV. The + states leaving the residual nuclei in their ground 
states are rather insensitive to the A currents for |q| < 500 MeV/c, since their excitation is mainly of 
longitudinal type. The 1 + states are instead excited via transverse responses and, even at 100 MeV, 
they are strongly affected by the presence of the A. The situation is intermediate for the 2 + states 
where both longitudinal and transverse excitations are important. 

The test of sensitivity of the cross section on the details of the correlation function has been done 
by using three different correlation functions taken from Fermi hypernetted chain calculations on finite 



nuclear systems [32, These are realistic correlation functions and differ only for few details. 

A first remark about our results is that the shape of the angular distributions of the cross sections 
remain essentially unchanged. This shape is insensitive to the correlations and it is determined by the 
decaying waves allowed by the angular momentum and parity composition related to the quantum 
number of the nuclear final state. We essentially confirm the arguments of Ref. [14] where a partial 
wave decomposition of the various nuclear final states is presented. 

We have already mentioned that the maxima appear to emission angles corresponding to the 
minimum values of the nuclear recoil momentum p r . The modification of the shapes due to the 
changes of the kinematics can be understood in terms of Fourier transform of the two-hole relative 
wave function with respect to p r . 

The SRC act on the size of the cross section. Calculations done with a set of schematic correlation 
functions, allowed us to understand the observed behavior in the results obtained with the realistic 
correlation functions. The size of the cross section is related to the function h{r) = f 2 (r) — 1 where 
/(r) is the two-body correlation function as it is usually defined. The cross section is more sensitive 
to the length of the healing distance, the distance where h(r) becomes zero, than to the value of h(r) 
for r = 0. An overshooting of the asymptotic value of f(r), corresponding to a change of sign of h(r), 
produces results similar to those obtained by reducing the length of the healing distance. 

The information about the SRC can be obtained only by a quantitative comparison between the- 
oretical predictions and experimental data. Qualitative features, such as the shape of the angular 
distributions, are not sensitive to the details of the SRC. Unfortunately a precise quantitative eval- 
uation of the (e,e'2p) cross sections is linked to the theoretical framework used to calculate it, and 
to the uncertainties in the required theoretical input. It appears clear that two-nucleon knock out 
experiments cannot be considered as the ultimate tool able to determine the characteristics of the 
SRC correlations. Instead, they have to be considered as another, useful and interesting, element of a 
puzzle, that together with elastic, inclusive and one-nucleon emission experiments should be described 
in a unique and coherent theoretical framework. 



APPENDIX A: Matrix elements of the one-body operators 



In this appendix we give the explicit expressions of the matrix elements ( pl|) needed to evaluate the 
response functions of Eqs. (|To|)-(13). Actually, what we calculate is the ^ 1 term as defined in Eq. 
( |37|) . The calculation of the various terms defining this transition matrix element proceeds by making 
a multipole expansion of the correlation function hij\ 



h 



'j 



H r ij) = h(n,rj, cos 9ij) 



EM 

L=0 



ri,rj 



Pl(cos 6i 



(42) 



where Pl are the associated Legendre polynomials. 



As shown in Eq. (32), the operator O n (q) is the Fourier transform of a one-body operator, in our 



case the electromagnetic operators (14) and fll5|). We perform a multipole expansion of the e~ ?qr 
term and, because we have chosen the z-axis parallel to the direction of the momentum transfer q we 
obtain 



Oo(q) = p(q) 



J 



J / d 3 rjj(qr)Y J0 (n)p(r) 



(43) 
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0±i(q) = J±i(q) = ^EH Jj 7 d3r 



J 



v7 



±i J (gr)Y|i )0 (fi) + (-i)yi J+ i(gr)Y^ il|; (n) 



where we have indicated with 



Yft,j(n) = E < A/^lr/lJM > Y XlM (n)e v 



J(r) , (44) 



(45) 



the vector spherical harmonics [^9| and with j\ the spherical Bessel functions. In the above equations 
we have defined I = \/2l + 1. 

From the above equations and the expressions (|i"4] ) and ([D]) we found convenient to express each 
multipole component of the one-body operator operator O v in terms of products of a factor depending 
on qr and another one depending only on the angular coordinates f2. For the charge operator (ry = 0) 
we have 

Oo(q) = E / d 3 rFj(qr) Oj^O) . (46) 



By taking into account Eqs. ( |i~4| ) and ([£]) we identify 



A r 



Fj(gr) = V^(-i) J Jjj(qr)J2 



G|(q)^3(!) + ^(q)^# 



and 



e>j, (ft) = y J0 (^) 



(47) 



(48) 



In the above equation we have indicated with Gg N (q) the electric form factor of the proton and 
neutron, respectively. 

For the magnetization current (?? = ±1) we have 



0±i(q) = E E F}(qr)Ol ±1 (n). 



(49) 



J " s=±l,0 

Taking into account Eqs. ( |l5|) and (]44| ) the factors in the previous equation are defined as: 



F s j(qr) = V2^(-i) J T(J,s)j J+s (qr) E 

i=l 

with G*]y[ N (q) the nucleon magnetic form factor and 

T(J,s) = 



G?(q) l + 73 (t) , G^(<J) 1-T3(* 
2rrii 2 



+ 



2m,- 



(50) 



-irjy/J + 6 s ,-i , s = -1 

J, 5 = 

irjyj , s = +1 



and 



(51) 



(52) 



oyn) = [y J+i (n)®^ . 

In the following, the coordinate where the one-body operators is acting will be labeled r\. The 
same notation has been used in Eq. (|37|). For the charge we define the integrals 



K(q;J;a,b)= J dnr? Fj(qn) <(n) R b (n) , 



(53) 



I(q;J,L;a,b,c,d) = J dnr? | dr 2 r 2 2 fc £ (n, r 2 ) Fj( 9 n) ^(n) # 6 *(r 2 ) i? c (n) R d (r 2 ) , (54) 
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H(L;a,b,c,<t) = J dnr? J dr 2 r 2 2 h L (r u r 2 ) R* a ( n ) R* b (r 2 ) R c (n) R d {r 2 ) , (55) 



and 



where R is the radial part of the single particle wave function defined in Eq. fl24j). Similar expressions 
are defined for the current operator by changing Fj(qr\) into Fj(qri) in the K and / integrals. 

The orthonormality of the single particle basis used to construct the Slater determinants | <3?f > and 
|3>(0, 0) > ensures that only the single particle wave functions directly involved by the electromagnetic 
operator and by the correlation function contribute to the matrix elements. 

Following the nomenclature of the diagrams shown in Fig. [3| we obtain 

4-7T 

(0)( 2 .i) = E I(q;J,L;p 1 ,p 2 ,h 1 ,h 2 )6 tpijthi 5t P2 ,t h2 ^ (—1) +v+ i pl ~ mpl+jp2 ~ mp2 

JLX L 



J L A W j Pl A j hl \ ( j P2 L j h2 

r] -N N-rj ) I -m pi -N + rj m hl J 1 -m P2 N m h 



(j P1 \\[Oj®Y L ] x \\j hl )(j P2 \\Y L \\j h2 ), (56) 

where 5 ay b is the Kronecker symbol. The expression of the (2.2) diagram is obtained from the above 
one by interchanging h\ and h 2 . The expressions of the (2.3) and (2.4) diagrams are obtained by 
interchanging pi with p 2 in the expressions of the (2.1) and (2.2) diagrams, respectively. 
For the three-point diagrams we have 

(0)3.i = EE K (T:J;Pi,i)H{L-i,p 2 ,hi,h 2 )5 tpi ,tJu,t hl ^t P2 ,t h2 

JL i 
Air 

-I_ (_l).?Pl- m Pl (_l)ip2 -m P2 (_l^ N +Ji- m i 

( jpi J ji \ ( j P 2 L 3h 2 \ ( 3i L jhi \ 
y -m pi 7] nii J y -m p2 N m h2 J y -rrn -N m hl J 

(j Pl \\Oj\\ji) (j P2 \\Y L \\j h2 ) (ji\\Y L \\j hl ) . (57) 

The expression of the diagram (3.2) is analogous to the above one but interchanging h\ with h 2 . The 
(3.3) and (3.4) diagrams are obtained by interchanging p\ and p 2 in the (3.1) and (3.2) diagrams 
respectively. For the other diagrams we have: 

(0)3.5 = EE K (l^ J ^h h 2)H(L;p 2 ,p 1 ,i,h 1 )5 Utth2 5 tp ^ t Jt P1 ,t hl 

JL i 
Air 

111 (_]V'i-m< (—\\N+j P2 -m P2 

L 2 

( ji J jh 2 \ ( j pi L j hl \ ( j P2 L ji \ 
y -nii rj m h2 J y -m pi N m hl J y -m P2 -N mi J 

(M\Oj\\j h2 ) (j P1 \\Y L \\j hl ) (j P2 \\Y L m ■ (58) 

Also in this case we obtain the expression of the (3.6) diagram from the above expression by inter- 
changing hi and h 2 . 

The above expressions are specialized for the charge operator. For the magnetization current, Oj 
must be changed into Oj^. The various charge and current transitions are calculated by substituting 



the above definitions in Eqs. (5€)-(|5q) and remembering that 



<./«Pi||./,,> = ' | J \ J aia + h + L), (59) 



<3a\\[Y L ®Yj] K \\j b> = (-1)L-J^=I L Q K q I <ja\\Y K \\jb >, (60) 
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<ja\\[Yj®*] J \\j b > = (-1) 



w a jajb'J I jb ja J \ tn l7 T \ 



(61) 



<ia||[Y Js ®<x] J ||j fe > = ("I) 



la+h+h+h Ml Xa + Xb + sJ + Ss,! 



\/47T + <W 




2 2 



+ lb + j + 1) , 



(62) 



where \j a >= \lb2ja >, with £(L) = 1 if L is even and zero otherwise and 



Xa = (-l) ia+> + ^(ia + i) = (la ~ ja)(2ja + 1) • 



(63) 
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